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Abstract 

Magnetosonic waves are intensively studied due to their impor- 
tance in space plasmas and also in fusion plasmas where they are used 
in particle acceleration and heating experiments. In the present paper 
we investigate the magnetosonic waves propagating in a multi-ion- 
species plasma perpendicular to an external magnetic field. Due to 
the presence of several ion species, this mode splits into two branches: 
high- and low- frequency modes. This opens a new channel of nonlinear 
interactions (between these two modes), and qualitatively changes the 
picture of turbulence in the long-wave region. Using the exact kinetic 
approach, we derive a general system describing the propagation of 
nonlinear ly coupled high- and low-frequency waves. This system in- 
cludes the KdV, Boussinesq, and Zakharov equations as limiting cases. 
Solitary solutions of the system of coupled equations are obtained. 
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I. DISPERSION RELATIONS FOR HIGH-FREQUENCY 
AND LOW-FREQUENCY WAVES 



We consider two-dimensional turbulence of magnetosonic (also referred 
to as magnetohydro dynamic or compressional Alfven) waves, propagating 
perpendicular to an external magnetic field in a coUisionless cold plasma with 
two ion species. As was shown in [1], in such a plasma the magnetosonic 
mode is split into two branches: low-frequency and high-frequency waves. 
The dispersion relation for them can be found from: 



and a marks the particle type: a — {a, 6, e} (we denote ion species as a 
and h) . This equation has two solutions, one of which corresponds to a high 
frequency wave and the other to a low frequency wave [1]. The dispersion 
law of the high frequency wave is: 
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and Va is a modified Alfven velocity: 
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For the low frequency wave Eq. (1) gives: 
where 



(7) 
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These relations were established in [1]. Dispersion relations (3), (7) are ob- 
tained in the long wave region and are valid when k « maxa=a,b{i^pa/ c} ■ 
We will also need the polarization vectors for each of these modes. In a sys- 
tem where the wave vector is directed along the a;-axis {ky = 0) we obtain: 



where 



e±{k) 
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,2 xI7jfe,-£±±) , 



(9) 



~g±{k) = g{i^±{k)) , 



(10) 



s±±{k) = e±{uj±{k)) 
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II. NONLINEAR CURRENT AND NONLINEAR RESPONSE TENSOR 



The nonlinear response Skkik^ is defined as: 

Skkik2 = ^*i{^)Sijiej{ki)ei{k2) , 



(12) 



and the tensor Siji can be found from the nonhnear current: 



jk,i = I Siji{kkik2)Ek,jEk^^i5{k - ki - ka) 

x5{u! — cui — UJ2) (Pki (Pk2 dcui duj2 ■ (13) 

In our hmit of cold plasma and small k this current can be obtained from 
the coUisionless kinetic equation, in which we neglect small terms kVTa/^ 
(where Vra is the thermal velocity of particles of type a) . A simple iteration 
procedure [4] then leads to the following expression for the nonlinear response 
tensor: 



^''^ ~ 47rm,(Q2_^2) {k2) - d^j) . (15) 

Particles of all types contribute to this tensor. As an example, let us calculate 
the nonlinear response for the high-frequency field Ei and low-frequency field 
E2. Consider the expression: 



k2s £^^\k) - Ssj ej{k2) 

in our system where k2y = 0. We will always assume that the field with 
subscript 2 corresponds to the low-frequency wave. Since for LJ2 = k2VA we 
have g « e±, we obtain that Gj{k2) = (0, —1), and: 



k2se^:;^e,{k2) = -ig^-h2 = -i^ . (16) 
Further, for high-frequency waves (a; = a;o) we have s± = g and therefore: 
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e*{k) = -^e^% 1), <h) = ^^-'Hh -1) , (17) 
where is the angle between the wave vector k and 

Now from (12), (14), (15) we get (with Bq being the external magnetic field): 
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Note that in a general case we should use an expression symmetrized with 
respect to the last two indices for the nonlinear response: 



Skk\k2 Skkik2 Skk2ki ■ 



Calculation of Skk2ki can be done in an analogous way: 



ICCUq 



e^^'t'-^^) j3k- 



(19) 



where 



We will suppose that (3ki « 1 and term (19) is negligible in comparison 
with (18). We will return to the meaning of this inequality in Section III. 

III. KINETIC DERIVATON OF THE NONLINEAR COUPLED EQUATIONS 

We use the kinetic approach for the derivation of nonlinear coupled equa- 
tions for the high- and low-frequency waves. Wc use Maxwell's equations for 
the fields and the coUisionless Boltzmann equation for the particle distribu- 
tion function. 

We represent the high-frequency field in the form 



where the field E"*" is concentrated at the frequency uj — uq, and the field E 
at the frequency lu = —luq. Maxwell's equations give: 



E = E+ + E 



(20) 




(21) 



where is the nonhnear current obtained from (13). This current is of 
the second order in field amplitudes, therefore if the field Ej in equation 
(21) is concentrated at cuq, the only term contributing to the current is: 
~ E^E^ , where E^ is the low-frequency field. In nonlinear current (21), one 
can also keep the higher order terms ~ E^E^E^ . If the amplitude of the 
low-frequency field is small, E^ ~ E^E^. one should keep the cubic term 
together with the quadratic one in (21). If the field Ej in (13) is a low- 
frequency one, the main contribution to the nonlinear current is given by 
E^E° and E+E'. 

Let us substitute Ej = ej{k)E^^k, where ej{k) is the the polarization 
vector of the field in (21), and multiply this expression by e*(k). Then we 
have for the high-frequency field: 



e*(k) \^-kikj + k%j - ^e,,j ej{k) = - — {e^-g) + -^ (22) 
and, expanding {e±_ — g) in small deviation of the frequency from ujq: 



£_L - fif ~ (a; - uo)-^{€± - g) 



{uj-uJo) , 



(23) 



we get: 



(a; - a;o - /.+fc^)^+, = -^J^e^j^ (k) . (24) 

This is the main equation for the high-frequency field. To close the equation, 
we express the nonlinear current in terms of the field amplitudes: 



jfik)^2JS,jkEnh)E^Mdl2 
+6JEijkiE+(k,)E+(k2)Er(k3) dl23 , (25) 

where we use the short-hand notation: 

dl2 = 5{u; — LOi — u;2)S(k. — ki — li2)duJidu;2dkid'k2 , 
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dl23 = 5{u! — LO\ — LOi — a;3)(5(k — ki — k2 — \i-^diO\dL02dLO-id^\6^i'ik.-}, . 

The response Si^k is symmetrized with respect to the last two indices, the 
the response T^ijki with respect to the last three indices. 

The second order response can be obtained using the general formula 
(13), which gives: 

'§(^^1^2) = -i^e'(^-'^^)(l + 0{i5k,)) . (26) 

We will consider the limit (3ki « 1 (i.e. /ciVa << -^^a) and will neglect 
the last term in (26). One can also show that the third order term can be 
safely neglected for our problem. 

Consider now the derivation of the equation for the low-frequency field. 
Multiplying equation (21) by the complex conjugate vector of the low-frequency 
field polarization, we get: 



= ^|^e:'»'i» . (27) 

This is the main equation for the low-frequency field. The contribution to 

the low-frcqucncy current is given by two terms proportional to E^E^ and 
to E^E^ . Consider the contribution of the first one. In this case, all three 
polarization vectors in formula (12) are related to the low-frcqucncy field. 
Using result (16) and the following expression for the polarization vector 
e(ki) = (sin 01, — cos0i), we get: 

qOO _ V ga(^Sin0i -iQaCOS(/)i) ^ .^ ^lak2 , . 




C.00 _ V gg ( -0; COS 01 - iOg sin 01 ) ^ ..'^^^ /oq^ 

'-'2 — 77^ K\ I, I'J-t: ■ v^yj 

47rma(si^ — UJ'^) >^la^2 

Multiplying by e*(k) = (sin0, — cos0), we get: 
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S'' = ^a-, ^^^^ - [u COS {4>-<l>,)+tn^sm{4>, -</>)] . (30) 

For cu « the second term in the brackets is small due to plasma 
quasineutrality, and we finally get: 



^ -^a-. ^ cos 0-01 . 31) 

Now let us find the high-frequency field contribution to the low-frequency 
nonlinear current. We choose the system where k2 = (A;2, 0) and get: 



e*(k) = (sin 0,- cos 0) , (32) 
e+(ki) = ^e-'^H^,-l) , 
e-(k2) = ^(-z,-l) . (33) 

The nonlinear response takes the form: 



S2^{k2k^) = -iS^+{k2ki) , 

and, multiplying by the polarization vector of the low-frequency field, we 
finally obtain: 



S''+-{kkik2) = ^(fcie'^'^^-'^) + /tae'^'^-^^^) . (34) 
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We are now ready to write the system of coupled nonlinear equations. We 
will consider the one- dimensional case. In this case we will use the follow- 
ing representation for the fields: E°(a;,A;) = {0, —l)E^{uj,k), E+(a;,A;) = 
■^{i,-l)E'^{uj,k), E~{uj,k) = -^{-i,-l)E~{uj,k). All of the polarization 
vectors are chosen to be fixed, i.e. independent of the sign of k. Intro- 
duce the envelope of the high-frequency field as follows: E^{uj,k) = 
E^{u, A;)e^*'^"*. Then, in Fourier representation this system takes the form: 



+ {third order terms) , (35) 



,,2 



^^)El, = -^JEiEtdl2 



1 OjJ- 



2Boc 

2B,vjn^.+^nE',dl2 . (36) 

We need to introduce some ordering to understand what terms can be ne- 
glected in (35-36). First of all, the most interesting case is when the disper- 
sion and the nonlinear terms in (35) are of the same order. This gives: 



- -^E^ , (37) 

where Vq is the characteristic velocity of the nonhnear wave. We also neglect 
the cubic term in (35). One can show that this can be justified only when 



-^\E I « -E . (38) 

Then in the left hand side of (36) we neglect the dispersion term ^^^fp^, where 

/i_ ~ This is justified, since this term is always small compared to the 

last nonlinear term in (36) due to (37). Thus we assume (37) and (38). If 
one then introduces, for convenience, the magnetic field perturbation Bi as 
E^ — Bi^, two self-consistent orderings are possible: 
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1. The first nonlinear term in (36) dominates the second one, i.e. Bi « 

E+- 

2. Both nonhnear terms in (36) are of the same order, i.e. Bi ~ E'^ . 
Consider the first ordering. Prom (37), (38) we get: 



Va 

Bo— » » Bi . (39) 
c 

Then, the hnear terms in (36) must be of the same order as the nonhnear 
terms. Comparing k'^E^ and the nonlinear term, we have: 

which is in agreement with (39). Consider the left half of inequality (39): 



f « ^ ^ 

Using (40), it is easy to check that this condition coincides with the previous 
assumption f3k << 1 (see (49)). Thus, in limit (39), we neglect the cubic 
term in equation (35) and the second nonlinear term in equation (36). 

Consider now the second ordering, Bi ~ In this case we keep both 
nonlinear terms in equation (36). Comparing these nonlinear terms with the 
linear ones, we find that they are of the same order if the velocity of the 
nonlinear wave Vq is close to Va- 

SV E+ 

— r.— ,dV=\VA-Vo\. 

VA -Do 

The condition pk « 1 together with (37) then gives: 



^ « S (42) 
Bo ' ^^^^ 

which is more restrictive than (41). We therefore demand (42), and the 
second ordering is then also self-consistent. In this case, we neglect the cubic 
term in (35) and keep both nonlinear terms in (36). 
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IV. SOLITONS IN A MAGNETIZED MULTI-ION PLASMA 



Consider the soliton solutions for system (35-36). Let us rewrite the 
system in the {x, t)-representation. For this purpose we introduce the "po- 
tential" $: = d^/dt. The system then takes the form: 



y dt dx^ ) cBq dx ^ ^ ^'^ 



As we have already seen, two limiting cases are possible, depending on the 
relative magnitudes of the nonlinear terms in (44). In the first case we neglect 
the second nonlinear term in (44) and get the following system: 





-— -^U = -^^|E+p (46) 

In the second case we keep both nonhnear terms in (44). In this case, the 
nonlinear wave has velocity Vq which is close to the velocity Va, and we can 
replace: ^ = in the right hand side of equation (44). Then, if we are 

interested in the wave moving in some particular direction (say, to the right) 
we can further simplify (44) by writing: 
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We then get the following system: 
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dt dx^ 



{vAdt^ dx) " AB^c ax ' ' + 2So • ^ 

Let us look for the stationary solutions of systems (45-46) and (47-48). 
It is convenient to represent the fields in the form: 



$ = $(^) , (49) 

where 6 = x — V^t, and Vq, A - some parameters. First consider system 
(45-46). Substituting expression (49) into it, we find: 



2^2 

EA + Eee^-^E^e , (50) 

C-Do 



Integrating equation (51) once, we get: 



where A is some constant of integration. Substituting this result into (50), 
we consider two cases: Vo < Va and Vq > Va- In the first case, the soliton 
solution has the form: 



^^vfWf^^' ^^^^ 
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where B = A — C = -^j^^^^^-^. The sohton solution exists only when 

B > 0. In the opposite limit Vq > Va, we have another soliton solution. 



E.J?^^^, (54) 

which exists only for S < 0. 

Note, that system (45-46) is analogous to the Zakharov system for coupled 
Langmuir and ion-acoustic waves [3]. The difference is in the sign of the 
nonlinear term in equation (45). Consider now system (47-48). Substituting 
(49) into (48) we obtain: 



where F is an arbitrary constant of integration. The soliton solution takes 
the form: 



_ ath{be) 
~ ch{hd) ' 



where 



6^ = -(A-Mi^)^-A , 

C-Do 

Va-Vq ^Dc 3c^ ~ 
Vo Bo M ' 



2 \%Blc^A^ 
a = 4 
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This solution depends on two arbitrary parameters, A and is valid for 
both Vq < Va and Vq > Va, is analytical at Vq = Va, and exists only for 
A < 0. 

V. CONCLUSION 

In conclusion, we have presented the general kinetic method for derivation 
of nonlinear equations for magnetosonic waves, propagating perpendicular 
to an external magnetic field in a multi-ion plasma. Equations (45)- (46) 
constitute the so-called Zakharov system [3] , first derived for Langmuir and 
ion-acoustic waves in non-isotermal (Tg ^ Ti) plasma. Without the high- 
frequency field and with dispersion term fi^k'^/VA, Eq. (44) is the Boussinesq 
equation for acoustic waves, and it can be reduced to the KdV equation 
exactly in the same way as we went from (44) to (48). 

Equations (45)- (46) can be used to describe weak turbulence of interact- 
ing high- and low-frequency waves. For this purpose, the equations should 
be rewritten in the Hamiltonian form, expanded in small amplitudes of in- 
teracting fields, and averaged over statistical realizations. Such a procedure 
is described in [5]. In [2] the same turbulence was considered using another 
method. The system can also be used to describe transition to strong turbu- 
lence. 

The system analogous to (47)- (48) was considered by Makhankov [6], and 
also by Nishikawa et al [7] in connection with near-sonic solitons of coupled 
Langmuir and ion-acoustic waves. 

We would like to stress that the developed method, though it is not as 
transparent as the hydrodynamical approach, is exact and rather straightfor- 
ward, and can be easily generalized to other plasma systems. Moreover, it, in 
principle, allows the consideration of such kinetic effects as Landau damping 
on a rigorous basis, which is impossible with the hydrodynamical approach. 

This work was supported by U.S.D.o.E Contract No. DE-AC02-76-CHO- 
3073. 
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